Helical spin-density wave in V^2-y02 



00 

u 
< 



X3 
O 

o 



> 

in 
o 

o 

oo 

ON 



o 



X 
S3 



T. Wolcnski, M. Grodzicki, J. Appel 
Universitat Hamburg, I. Institut fiir Theoretische Physik 
Jungiusstrafie 9, D-20355 Hamburg. 
(Submitted: November 6, 1997) 

Recent neutron scattering and nuclear magnetic resonance experiments have revealed that the low 
temperature phase of doped V2-j;03 is an itinerant antiferromagnet with a helical spin structure. 
We use a band structure calculation as the point of departure to show that these experiments are 
in agreement with mean field results for an Overhauser spin-density wave state. The influences of a 
finite life-time and of dilute magnetic impurities are discussed. 



I. INTRODUCTION 

Stoichiometric V2O3 is generally considered a classic 
Mott-Hubbard system displaying a first order transition 
at Tat = 155 K from a paramagnetic metal (PM) to a local- 
moment antiferromagnetic insulator (AFI). It is com- 
mon wisdom that this transition is driven by the strong 
Coulomb interactions between conduction band electrons 
which opens up the Mott-Hubbard gap. The antifer- 
romagnetism, which lifts the degeneracy in the insulat- 
ing phase, is therefore a consequence of the interaction- 
driven opening of the gap between the upper and lower 
Hubbard band. 

Experimentally the situation appears to be quite differ- 
ent in doped V2-j,03, where the low temperature phase 
is an antiferromagnetic metal (AFM) with an incom- 
mensurate Overhauser spin-density wave (SDW), similaj: 
to metallic chromium. Neutron scattering experimentsu 
show that the SDW has a nesting vector Q = 1.70 c* (c* 
is a reciprocal lattice vector in the hexagonal unit cell of 
V2O3), and that the transition appears at a T/v of about 
10 K; both T/v and Q are almost independent of dop- 
ing in the regime probed. The ordered moment is only 
^ = 0.14: as compared to the 1.2 in the 

Moreover, the magnetic excitations are founda to have 
a bandwidth of about 20 fc^T/v. In a local moment pic- 
ture one would expect a magnetic bandwidth of order 
J ^ UbTn. The large bandwidth can be understood 
in the framework of SDW theory, where, however, ex- 
citations are found on electronic energy scales. An in- 
teresting observation is that even close to the metal- 
insulator transition the magnetic excitations correspond 
to the SDW spin structure rather than the magnetic or- 
der in the antiferromagnetic Mott insulator phase. This 
suggestsu that the local vanadium moment formation in 
the AFI phase is accompanied by a change in the elec- 
tronic configuration. l-. 

NMR experiments by Langenbuch and Piepeiu have 
confirmed the neutron scattering results for the ordered 
moment of about 0.14 /is and the orientation of the SDW 
along the c-direction. The magnitude of the local mag- 
netic field is identical on all V-sites, compatible with a 
helical rather than linear polarization. At lower doping 
the temperature dependence of the spin-echo intensity 
also shows T/v —10 K. However, at doping y > 0.038 the 
transition observed with NMR shows T/v — 50 K, with 



no special feature in the V relaxation rate around 10 K. 
The phase diagram combining these experimental results 
is shown in Fig. An additional important contribution 
of NMR is the study of the minority V sites neighboring 
the vanadium vacancies in c-direction. These sites are ap- 
proximately in a V*+ (S=l) state. Furthermore, the 10 K 
transition could be intrinsic to the V2O3 system, while 
the magnetic impurities increase the transition temper- 
ature above a certain impurity concentration. We will 
discuss this scenario in Sec. 0. The homogeneous suscep- 
tibility at higher dopingo is very sensitive to whether the 
measurement is made in zero field or in the presence of 
the earth magnetic field. In the former case, for ?/=0.052 
a transition is seen at 56 K, while there is no feature at 
10 K. In the latter case, there is only some hysteresis left 
at higher temperatures, while there is a pronounced peak 
at 10 K. 

The order of the phase transition is an issue not com- 
pletely settled by these experiments. While the neutron 
scattering data are consistent with a BCS-like second or- 
der transition at 10 K modified by fluctuations, the NMR 
data show, at least for the 50 K transition, a large dis- 
continuity in the local magnetic field at T/v. 

There have been attempts to qualitatively study the 
itinerant AFM phase using the LISA (local impurity self- 
consistent approximation) method starting from the limit 
of strong correlationgj, where the key factor for the clos- 
ing of the Hubbard gap is the magnetic frustration caused 
by the defect sites and the V^"*" moments. However, as 
pointed out above, the experiments indicate a conven- 
tional Overhauser SDW resulting from Coulomb interac- 
tions at nested Fermi surface sheets. In this paper we 
show that it is possible to explain the experiments with 
a molecular-field Hartree-Fock SDW calculation starting 
from accurate band structure results and taking into ac- 
count realistic electron-electron scattering and multiple 
bands at the Fermi level. 

In the following section we discuss the band struc- 
ture calculation, which is the point of departure for the 
Hartree-Fock SDW calculation with life-time effects dis- 



cussed in Sec. HI. The results will be discussed in Sec. IV 



In Sec.0 we compare our results with the NMR and neu 
tron scattering experiments and discuss the importance 
of magnetic impurities. 
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II. BAND STRUCTURE 

The first band structure^ calculations for V2O3 date 
back to the early seventies.El However, at that time only 
a tight-binding calculation with a rather Kmall set of or- 
bitals was feasible. Recently, MattheisaO performed a 
state-of-the-art calculation using the scalar-relativistic 
linear augmented-plane-wave (LAPW) method. Wf 
make use of his results in the form of a Slater-Koste: 
calculation using the parameters given by Mattheiss. 

In the Slater-Koster method the band structure is con- 
structed from parameters obtained at special points of 
high symmetry using LDA-APW results. Then, the 
energy eigenvalues at general k are interpolated using 
a tight-binding-like formalism, thus preserving the full 
symmetry of the lattice. 

To find the energy eigenvalues, in general one has to 
calculate integrals of the form 



J 0;(r-R^)i70^,(r- (1) 



for a Hamiltonian 



(2) 



If ^, ^' and K are all different, the integrals are three- 
center integrals; if k = /i or k = /i', they are two-center 
integrals. The three-center integrals tend to be smaller, 
and the angular dependence differs only for the I > 
contribution, so that it is a good approximation to only 
consider two-center integrals (extended Hiickel approx- 
imation). The wave functions and the potential 
are not very well known. Therefore, the concept of the 
Slater-Koster method is to treat the two-center integrals 
as parameters obtained from LDA calculations at special 
points, rather than calculate them microscopically. The 
remaining non-trivial task is to construct the full matrix 
elements between all 44 orbitals considered in the trig- 
onal unit cell of V2O3: Five 3d-orbitals for each of the 
four vanadium atoms and one 2s- and three 2p-orbitals 
for each of the six oxygen atoms in the unit cell. 

In Fig. H we show the bands in a trigonal crystal struc- 
ture plotted along lines of high symmetry. The sym- 
metry points are labeled using the hexagonal notation, 
the Fermi energy is set to be zero. In Fig. || we plot the 
density of states in the vicinity of the Fermi energy, stem- 
ming predominantly from the V(3d)-bands. We also plot 
the density of states from the two bands (labeled band 
A and B respectively) which contribute about 80 percent 
of the total density of states at the Fermi energy. Note 
that the densities of states of band A and band B at the 
Fermi energy are almost equal. 

The Fermi surface is rather complicated, and there- 
fore we restrict ourselfs here to showing two cuts through 
the Brillouin zone: Fig. |4| a) is in the a-b {kz=0) plane 
and reflects nicely the trigonal symmetry of the lattice: 
a three-fold rotation axis in c-direction combined with 
inversion symmetry around the origin. This symme- 
try leads to magnetic frustration in the a-b plane, in 



agreement with the observed SDW nesting vector in c- 
direction. In Fig. ^ b) we show a second cut in the a-c 
{ky=0) plane, where we indicated the main nesting fea- 
ture of the Fermi surface with an arrow. The nesting vec- 
tor has been calculated numerically from the maximum 
in the susceptibility (cf. Sec. IV), the arrow in this figure 
only serves for illustration purposes, in particular since 
the Fermi surface is truly three-dimensional and thus a 
single two-dimensional cut can be misleading. The nest- 
ing part of the Fermi surface we indicated extends into 
the 6-direction, thus giving a sufficient density of states 
to produce a spin-density wave as will be shown in the 
next section. 

Strictly, the band structure calculation is performed for 
stoichiometric V2O3. However, the doping in the sam- 
ples experimentally probed is sufficiently small (y < 0.05) 
that even at the highest doping levels the chemical po- 
tential is only shifted by a few meV, justifying a rigid 
band approximation. 



III. HELICAL SPIN-DENSITY WAVE THEORY 
A. Molecular-field gap equation 

A helical spin-density wave instability, where electrons 
with momentum k and holes with momentum k+Q and 
opposite spins pair below a critical temperature and thus 
produce a gap at parts ofihe Fermi surface, has first been 
proposed by OverhauserQ for the Hartree-Fock electron 
gas. While for the three-dimensional electron gas this 
instability is prevented by screening, certain Fermi sur- 
face topologies with sizeable parallel parts of the Fermi 
surface connected by a single vector Q (nesting) lead to 
the appearance of an SDW phase in actual materials like 
chromium£3 

Starting from a Hubbard Hamiltonian 

-ff = ^ ekCk_„Ck,Q + J7 ^ c|,.^Ck',a'Ck'-q,a'Ck+q,a, 



k,k',q a, a' 



(3) 



the assumption of a finite expectation value for Sq = 

'^ki'^+Q^'^U'^^^'^ ^^'^ Hartree-Fock approximation 
leads to a molecular-field equation for the SDW gap by 
means of a standard Boguliubov transformation, rather 
similar to the BCS gap equation. Since the nesting vec- 
tor Q is determined by the Fermi surface topology, it will 
in general be incommensurate with the lattice. The gap 
equation has beeii-generalized to the incommensurate sit- 
uation by Youngjli in our notation it reads 
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^ V(fk-ek+Q)2 + 4A2 
The quasi-particle energies are given by 

i^k" - ^^^^ ± 



(4) 



(5) 
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and /(£■) = l/(exp(i;/fcsT) + 1) is the Fermi function. 

The assumption of a finite expectation value Sq leads 
to a spiral polarization of the SDW,lii as we can show 
from considering Sq = Sq = Sq and SZq = S^q = 
Sq. As usual, they are related to S^ and 5^ by S^ = 
{S+ + S-)/2 and Sy = {S+ - S-)/2. The Fourier trans- 
forms then are 

5^(R)-(5+(R) + 5lp(R))/2 

= 5Q(exp(iQ • R) + exp(-iQ • R))/2 
5q cos Q • R (6) 
5^(R) = 5QsinQ-R, (7) 

which is a spiral polarization. If the nesting instability 
were strong enough to simultaneously support an SDW 
with opposite helicity, we would get a linear polariza- 
tion. In the present case, however, the subtle modifica- 
tion of the Fermi surface through the formation of one 
spiral SDW is sufficient to prevent this additional, oppo- 
site SDW. 

B. Life-time effects due to quasiparticle scattering 

The analysis we presented so far is for a one-band 
model without any scattering. As we show in the pre- 
ceding section, at least two bands are relevant for the 
low energy physics of V2O3. Testing the nesting proper- 
ties of these two bands, it shows however, that only one 
of them shows significant nesting leading to a possible 
SDW instability. Hence, our approach is to assume that 
the SDW results from the nesting band A, while band B 
enters the physics via scattering processes from the part 
of the inter-band Coulomb interaction not already taken 
into account in the LDA band structure. Since the den- 
sities of states at the Fermi energy of bands A and B are 
comparable, this scattering is expected to significantly 
shorten the life-time of band A quasi-particles. 

The scattering results from both electron-electron in- 
teractions, TpT^, and the scattering of electrons from the 
vanadium vacancy sites, t^^. The sizeable electron- 
electron scattering, which has the formll3 (Na and Nb 
are the densities of states for the respective bands) 



cNA{0)NB{0)iTr\kBTf + {hcof), 



(8) 



is responsible for the term in the electrical 

resistivity,tj while Tj^p determines the residual resistivity 
at T = 0. We thus have 

r-i = cNAiO)NBm^'ikBT)^ + (f^^f) + ^i7„p- (9) 

The constant c contains the scattering cross section. To 
estimate c we use the relation between resistivity and 
scattering rate. 



P 



(10) 



The carrier density n can be estimated from the Hall 
coefficient. 
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(11) 



Experimentally,^! Rjf ^3.5x10 cm"^ C ^ , resulting in 
a carrier density n = 1 .8 >jJ|0^^cm'^ . The coefficient of the 
term in the resistivityy is 6.0 x 10"^ cm K"^. Tak- 
ing all these values derived from different experiments at 
face value, a realistic effective mass m* of a few times m 
results in a scattering rate hT~^ /Ub ~ 10 K at T/y. One 
should consider this an estimate of r"^, which serves to 
prove the consistency of our approach, rather than an ex- 
act value, and hence we will parameterize the scattering 
rate by multiplying with a scaling parameter a, which 
varies between 0.1 and 1, i.e. 



Tit = a{kBTy 



(12) 



To take these life-time effects into account we make 
use of an approach first applied to the case of BCS 
superconductors with strong damping of the quasi- 
particles by Bardasis and Schrieffer.tS These authors 
derive the molecular-field gap equation including the 
Nambu-Eliashberg function Z(q, w), which is related to 
the scattering rate via 



(13) 



Z(k, Lu) = ReZ(k, lu) [ 1 + i 



The real part of Z(q, uj) renormalizes the band structure 
and is already taken into account in LDA. Thus, we set 
RcZ(k, cj) = 1. The imaginary part of Z(q, lj) is taken 
into account by straightforward generalization of the gap 
equation of Ref. |l^ to the case of an incommensurate 
SDW and to finite temperatures by employing the usual 
summation technique over Matsubara frequencies; one 
arrives at 



N ^ 



^V(fk-ek+Q)2+4A2 



dio Z 



(14) 



IV. RESULTS 

The necessary first step is to determine the nesting 
vector Q. There are two basically equivalent ways to do 
this: One can examine the q-dependent susceptibility 



X°(q) = E 



/(£k+q) - /(£k) 
fk+q — fik 



(15) 



and look for a maximum at a position qmax- Within 
RPA, where the full susceptibility is given by XRPA(q) = 
X°(q)/(1 ~ ^^X°(q))7 with increasing U the phase transi- 
tion will first occur for Q = Qajax- Equivalently, one can 
study the gap equation Eq. (|lj) , and calculate the value 
of q that leads to the minimum value of U necessary to 
satisfy the gap equation at a given temperature, i.e. to 
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produce the SDW phase transition. Both methods lead 
to the same result Q = (1.69 ± 0.03) c* in surprisingly 
good agreement with the neutron scattering results. EJ In 
Fig. H we show the susceptibility plotted for different val- 
ues of q along the c-direction in the paramagnetic state 
close to the SDW transition. 

Since the effective screened interaction U is not known 
very well from first principles, we use it sis-^ parameter 
to fix Tm- From optical conductivity dataO U has been 
estimated to be about 1 eV, a value also used in the 
LISA calculationaJ. In Fig. || we plot the interaction Uc 
necessary to obtain a given critical temperature for the 
SDW transition. 

The next step is to solve the gap equation for given 
U and r for the SDW gap A(r). In Fig. | we show the 
results for a of 10 K and different scattering rates. 
Clearly, the jump in the order parameter at T/v increases 
with increasing scattering rate. However, even at the 
lower limit of the reasonable range of scattering rates 
this discontinuity is sienificant, in agreement with the 
findings of Pieper et alu. 

There has been some discussion in the context of 
chromium as to how a mean-field SDW transition can be- 
come first o*der. One suggestion is the effect of a "reser- 
voir band" jl3 which is polarized by the formation of the 
SDW. Clearly, a strong coupling to the lattice will also 
enhance the tendency towards a first order transition. In 
the present case, however, the strong lifetime effects of 
our two-band model already produce a first order transi- 
tion. The life-time effects are stronger for the 50 K tran- 
sition than for the 10 K transition, so that it is possible 
that the 10 K transition is modified once fluctuations are 
taken into account. 

It is hard to estimate the size of the ordered magnetic 
moment from this calculation, since we did not test the 
full q-dependence of the order parameter. However, we 
calculated that about ten percent of the Fermi surface 
make the main contribution to the SDW instability, in 
agreement with an estimate of Bao et aln. 



energy E-mt = Hj • Sj, where Hj is the local field due 
to the SDW and Sj the impurity spin. If the RKKY- 
type interactions between the dilute impurities are strong 
enough, this mechanism will also lead to the formation 
of a domain structure, which could explain the difference 
between neutron scattering and the more local NMR ex- 
periments. 

An open question is whether the delicate field depen- 
dence in the experiments of Pieper et aln is related to 
an ordering of the V*^ moments, preventing the forma- 
tion of helical SDW domains. Problematic is that the 
critical field, which would have to overcome the RKKY 
interaction between the dilute moments, is only of the 
order of 1 G. It is known-Lhat strong non-linearities ex- 
ist in spin glass systems,Lj where comparable fields are 
sufficient to suppress a sharp cusp in the susceptibility. 
However, this point clearly requires further attention, to 
clarify the experimental findings. 

Another open question is how the Mott transition at 
155 K, which apparently is accompanied by a change in 
the electronic configuration — and also in the lattice sym- 
metry, which becomes monoclinic in the AFI insulator — 
is influenced by taking into account the realistic multi- 
band bandstructure. Particularly the intriguing experi- 
mental observation, that the fluctuations on the metallic 
side of the metal-insulator type are of the SDW rather 
than local moment type with a magnetic structure given 
by Fermi surface topology, show that for the physics of 
the Mott transition in V2O3 both strong correlation and 
electronic structure effects are important. This problem 
is under investigation. 
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V. DISCUSSION AND CONCLUSIONS 

Taking an LDA band structure calculation as point of 
departure, we find a SDW transition at 10 K with re- 
alistic parameters with a nesting vector q=1.69 c* and 
a helical spin jStructure in agreement with neutron scat- 
tering results.El We believe that this SDW is the intrinsic 
low temperature phase of doped ¥2-^03. This leaves the 
puzzle of the partial disagreement of neutron scattering 
and NMRq as to what the transition temperature of the 
higher doped compounds is. 

We suggest the following scenario for the higher doped 
compounds. The V site neighborinE-.a vanadium va- 
cancy in c-direction is in a V^^ stateo with n = I^b- 
These moments are very dilute, at the relevant concept 
trations y < .05. Overhausenlj and recently WernerE2l 
have shown how an incipient SDW combined with mag- 
netic impurities can lead to an SDW transition. The 
SDW is driven by the maximization of the interaction 
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FIG. 2. Bandstructure in the vicinity of the Fermi energy 
(e = 0) plotted along lines of high symmetry of the hexago- 
nal-type Brillouin zone. 




doping y e [eV] 



FIG. 1. Phase diagram of V2-y03 as derived from neu- 
tron scattering, susceptibility, and NMR measurements. PM: 
paramagnetic metal, AFI: antiferromagnetic insulator, AM: 
antiferromagnetic metal, PI: paramagnetic insulator. 



FIG. 3. Total density of states near the Fermi energy and 
partial density of states for the two bands, that contribute 
most at the Fermi energy. 
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FIG. 4. Fermi surface in the Brillouin zone of V2O3. a) 
Plot in the fe^j-fcy-plane (fe^^O). b) fc^j-fc^-plane {ky=Q). The 
arrows indicate regions of parallel Fermi surface sheets of band 
A connected by the numerically determined nesting vector Q. 
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FIG. 5. The RPA susceptibility for U = 0.4 eV without any 
quasi-particlc scattering in the paramagnetic state just above 
the SDW transition. Other values of q/c are accessed using 
the periodicity and the inversion symmetry of the lattice. 
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FIG. 6. Critical interaction necessary to achieve a given 
Tn- The electron-electron parameter a determining the scat- 
tering strength is explained in the text. 
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FIG. 7. Order parameter A(T) plotted against reduced 
temperature T/Tn for different strengths of quasi-particle 
scattering. 
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